Frequency dispersion of photon-assisted shot noise in mesoscopic conductors 
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We calculate the low-frequency current noise for AC biased mesoscopic chaotic cavities and diffu- 
sive wires. Contrary to what happens for the admittance, the frequency dispersion is not dominated 
by the electric response time (the "RC" time of the circuit), but by the time that electrons need 
to diffuse through the structure (dwell time or diffusion time). Frequency dispersion of noise stems 
from fluctuations of the Fermi distribution function that preserve charge neutrality. Our predictions 
can be verified with present experimental technology. 
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Charge neutrality in mesoscopic devices is enforced by 
Coulomb interaction on a time scale {trc) given by the 
product of the resistance and the capacitance scales. The 
typical time (tjj) for non interacting electron motion is 
instead fixed by either the dwell or diffusion time, de- 
pending on the transparency of the interfaces between 
the system and the electrodes. The typical physically 
realized situation is Tf/c td- electrons, which would 
normally slowly diffuse, are pushed to run by the electric 
fields that they are generating themselves by piling up 
charge. The consequence is that the typical response time 
of the device is = + t^^ '^rc- "^^^^ been 
shown for the frequency dependence of both the admit- 
tance 0,0 and the noise [1,0, [Ml^ in mesoscopic chaotic 
cavities and diffusive wires. The inverse of the diffusion 
time appears instead as the relevant energy scale for the 
voltage or temperature dependence for both the conduc- 
tance and the noise in superconducting/ normal metals 
hybrid systems 0jBS^3- Indeed, in this case the en- 
ergy dependence is due to interference of electronic waves 
that does not induce charge accumulation in the systems. 
To our knowledge, in normal metallic structures a disper- 
sion on the inverse diffusion time scale has been predicted 
so far only for the third moment of current fluctuations 
[3 and for the finite frequency thermal noise response 
to an oscillating heating power An alternative, and 
less investigated possibility, is to study the low frequency 
current noise as a function of the frequency f2 of an exter- 
nal AC bias. The noise for a quantum point contact was 
calculated a decade ago and later measured 
Since the quantum point contact is very short, electron 
diffusion does not introduce any additional time scale 
in the problem and the resulting frequency dispersion is 
simply linear More recently, the noise for an AC 

biased chaotic cavity was considered [isl llij in the limit 
of small fields eV/TiQ, <C 1 (V is the amplitude of the 
AC bias and e is the electron change). The authors of 
Ref. ^3 found that the noise disperse only on the r scale, 
that is the combination of the diffusion and the electric 




FIG. f : Schematic of the system. A chaotic cavity connected 
to the electrodes through arbitrary coherent connectors char- 
acterized by a set of transparencies {Tn}- 



response time. This should be contrasted with the result 
of Ref. |l7| where the noise in a diffusive wire was studied 
when the conductor was biased by a short voltage pulse. 
Even if the dispersion of the noise is not considered in 
that paper, the authors find that the ratio of the diffu- 
sion time to the time duration of the pulse may affect 
the observed noise. This indicate that a dispersion of the 
noise on the external frequency could be present. 

In this paper we derive analytical expressions for the 
photon-assisted noise in chaotic cavities [Eq. ^V6\ in the 
following] and numerical results for diffusive wires. We 
show that the noise depends on the AC frequency f2 
on the diffusion time scale. We take into account the 
Coulomb response that guarantees charge neutrality over 
times longer than trc- Within our calculation, the ad- 
mittance has no structure at « 1 /td , while the deriva- 
tive of the noise with respect to the frequency shows a 
clear maximum. 

We start by considering a chaotic cavity (see Fig. ^ 
connected to two metallic leads through two barriers 
characterized by a set of transparencies {T^} where k 
takes the value 1 or 2. We define the conductances Gk = 
Gq the total conductance G = GiG2/(Gi + G2), 

and a dwell time to = 2TT'hGQ/{Gi + G2)S, where S is 
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the level spacing of the cavity and Gq = e^/(27r?i). We 
assume that the dwell time is much shorter that the in- 
elastic time. We also assume that = Gi + G2 ^ Gq 
so that we can use scmiclassical theory to describe trans- 
port and neglect Coulomb blockade effects. The cavity 
is coupled capacitively to the leads through the capaci- 
tances Cfc and to a gate through Cg. One can then define 
the typical electric response time as tuc = / G where 
Ce = Ci + C2 + Cg. If linear dimension of the cavity L 
are much larger than the Fermi wavelength Xp, the ratio 
trc/td « S/Ec « (Xf/L)''-^ < 1 where Ec = eVCs 
is the Coulomb energy and c? = 2 or 3 is the dimension 
of the cavity. Three different time dependent voltage bi- 
ases are applied to the gate and the two contacts {Vg, Vi, 
and V2), clearly the current depends only on two voltage 
differences, we keep the three voltages to simplify the no- 
tation. The cavity resistance is negligible with respect to 
the contact resistances so that we can assume a uniform 
electric potential inside. 

We begin by considering a very simple model of charge 
transport, where diffusion and electric drift are described 
classically. The charge Q in the cavity is related to the 
electric potential V of the cavity itself by the relation 
Q = J2k Ck{V~Vk), where fc = 1, 2, and g. The equation 
of motion for Q is then 



dtQ + Q/TD + Ii + l2 = 0, 



(1) 



where Ik = Gk{V — Vk). The equation can be conve- 
niently solved in terms of V{uj), the Fourier transform of 
the cavity potential: 
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where we defined ak = Gk/G^ for k 
When Tjic ^ tjj, the response time r = trctd / (trc + 
Tjj) is simply given by tuc- The admittance has a similar 
frequency dependence. For trc ^ 'Td and cjt <C 1 one 
thus finds the very simple result 



(2) 



and Ii{uj) = —l2{uj) = GT,aia2{V — Vi): charge neutral- 
ity is perfectly enforced for low frequency drive. As it 
follows from Eq. (Q, only Ik is finite, since Q and dtQ 
vanish in this limit. In the following we will consider 
this experimentally relevant low driving-frequency limit 
^trc ^ Ij which will enable us to use the uj <^ ^/tuc 
response for V{uj). 

Let us now consider the noise. We use the scmiclassical 
description of non-equilibrium tra nsport p rovided by the 
time dependent Usadel equations 
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DViGVg) + [ihidt, +dt,) + eV{h) - eVit2)] 5 = 0, 



where D is the diffusion coefficient in the cavity and the 
Keldysh Green's function satisfies the normalization con- 
dition: 
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5^,5{h~t2). (4) 



Since the conductance is controlled by the contacts, 
G{ti,t2) is uniform in the cavity. Eq. Q can then be 
integrated on the volume of the cavity, and relate Q to 
the current through the interfaces . The equation for 
Q takes then the form of a charge conservation equation: 

Y^ik- iTD [Mdt, + dt,) + eVih) - eVit2)] G ^ 0, 

k 

(5) 

where the voltage in the cavity is given by the electroneu- 
trality condition (O and 
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are the spectral currents that relate the Green's function 
of the electrodes Qk to the Green function of the cavity 
|2ll | . The form of Gk is given by the equilibrium solution 
for metallic leads of Eq. explicitly, in presence of a 
time dependent potential, one has: 



Gok{ti,t2) 



5{h-t2) 2Fk{tiM) 
-5{ti-t2) 



(7) 



where hBt^t) = el4W, Fk{tiM) = e'*''^''^ F,g{ti - 
t2)e-*'^'=(*2), and i^egW = /(d£/27r)e*^*/'"nanh(e/2T)] = 
~iVT / sinh(7rri/?i), with T being the temperature 
{ks = 1). 

In order to calculate the noise we need to modify the 
Green's function of the electrodes by introducing two 
counting fields |2(1.E3|: 



Gk = e'''-^''/^Goke 



(8) 



(3) 



where Ux is a Pauli matrix. Solving Eq. Q together 
with Eq. (O gives the value of G from which, using 
Eq. lO one can derive the counting field dependent cur- 
rent: 4(i,Xi,X2) = GkTY{axik{t,t)}/2. Derivatives 
with respect to xi ~ X2 give all the moments of current 
fiuctuations. We assume that the potential difference 
between the two leads is harmonically oscillating at fre- 
quency VL. Then the time average over the period 2tt/Q. 
of Ik{t) and ~'idlk{t)/d{xi — X2) gives the DC current 
and the low frequency noise, respectively. 

Let us first discuss briefly the current. For vanishing 
Xfe, G has the form given by Eq. ^ and it depends on 
a single function: F(ti,t2). It is convenient to change 
gauge by defining F{ti,t2) = e"*"^*^^ F{ti,t2) e-"''^*-'^ 
where fidt4>{t) = eV{t), and introducing new time vari- 
ables = t2 ~ h and t+ = (ti -I- t2)/2. F{t+,t^) is 
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a periodic function of t+: F{t+,t_) = e™"*+F„(t_) 
and its explicit expression in terms of the Fourier com- 
ponent of F'^'^^ comes from the solution of the equation 
of motion 0: 



Fn{t- 
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(9) 



Here Fi''\t) = Feq{t)Jn {2AkSiii{Qt/2)), J„ are Bessel 
fimctions, Ak = e{V^ - V°)/{hn), and V° are the am- 
plitudes of the AC field: Vfe(i) = Vfe°sin(m). In order to 
make a connection with the classical description given by 
Eq. we relate F to the charge on the cavity: 



Q{t+)^^{e/S)F{t^^O,t+). 



(10) 



Usadel equations lO for ^ reduces then to the 
continuity equation ^ for the total charge. Note that 
when the voltage time dependence Q is enforced, Q van- 
ishes identically, as it can be verified by calculating the 
limit t_ ^ in Eq. On the other hand, the energy 
distribution function, F(t+,e) = J t+)e*^*-/'^, 
varies periodically in time, and its dependence on the AC 
driving frequency, O, is on the scale of the inverse diffu- 
sion time T^^. When an electron enters the cavity, after 
a very short time ~ trcj the charge rearranges to keep 
the cavity neutral, the distribution function, instead, will 
relax on a much longer time given by t^. Since zero- 
frequency noise probes the electronic distribution func- 
tion, it should depend on the frequency on the same 
scale. 

To obtain the noise, it is sufficient to expand the 
Green's function up to first order in x: G = Go + Gii{xi ~ 

a )/2 -I- . . ., and substitute this expression into Eq. (0) 
llfl l . The term Go is given by Eq. ((JJ) with F given by 
Eq. lO. Note that in principle also V{uj) could depend 
on Xk- By going back to the Keldysh action formulation 
we verified that the corrections are of order higher than 
one in Xk, they thus do not contribute to the calculation 
of the noise. In order to fulfil the normalization condi- 
tion El we use the parametrization for Gi proposed in 
Ref. Em . 
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By defining X ^ Xi - X2 and choosing Xi = a2X, X2 = 
—aix we obtain W' = and the kinetic equation for W: 

iil + TDdt^)W = ^Xfe«fc[(l-/3fc)(^°^ + ^fe°^fe) 



f3kiF o Fk + Fk o F)] , 



where f3k = Gg/Gk E„ ^^''(1 - ^i^'^) are the Fano fac- 
tors of the two junctions and with o we indicate time 
convolution. To obtain the low frequency noise we need 



(12) 




FIG. 2: Frequency dependence of a differential photon- 
assisted shot noise in the symmetric chaotic cavity (Gi = G2) 
under fixed flux A = e(Vi° - V2)/{hn). The magnitude of 
dS/dO. is normalized to its value at f2 = 0. Curve (1) A = 1.0, 
(2) A = 2.0, (3) A = 3.0, (4) A = 4.0, (5) A = 5.0. For sym- 
metric cavity the curves appears to be independent of the 
transmission distribution of the contacts. 



W averaged over one period, the n = Fourier compo- 
nent of Eq. ifT^ thus suffices for our purpose. Substitut- 
ing the expression for W into the current definition and 
evaluating the traces we obtain: 



S = G^Y. 



hujoil + r) coth 



"o(i+r) 
2T 



1 + uolr^n^ 



[TJ{A^,A2) + J(Ai, Ai) + ^2 J{A2,A2)] 
+ 2Gsr(l-/3ia2-/?2«i), (13) 

where J{Ai,A2) = J„+i(Ai) Jz(Ai) J^„„(-A2) Jr(-A2), 
Ti = (ai + Pia2 - T)/2, T2 = (^2 + Piai - T)/2, T ^ 
ot\a2 + /Siai + /32Cki, and the amplitudes are A\ = a\A^ 
and A2 = -a2A, with A = {V° - V2)/hn. Expression 
(|13|l is one of the main results of this paper. For flru ^ 1 
it reduces to Lesovik-Levitov result [13I with the effective 
Fano factor, J-", appearing at the place of the quantum 
point contact Fano factor. For small A we recover the 
result of Ref. 0| : at A^ order there is no dispersion. For 
all other cases a frequency dispersion is present, as it can 
be seen from Fig. 13 where dS/dft is shown as a function 
of for small temperature T hn. In particular, we 
find that dS/d£L displays a weak maximum for f2 ~ I/td 
reminiscent of the reentrant behavior in superconductors. 

The driving frequency dependence of the photon- 
assisted noise may also be generated by electron diffu- 
sion in disordered metals. To illustrate this point we 
consider a one-dimensional diffusive wire of length L 
and dimensionless conductance g = G/Gq ^ 1. We 
employ the same procedure used for a chaotic cavity, 
with the important difference that we have to take into 
account the spatial dependence of G in the wire de- 
scribed by Eq. ©. The condition jSjl is substituted by 
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the neutrality condition on the potential along the wire: 
eViz,t) = eV°sm{nt){L - z)/L, with < z < L. We 
assume now that contact resistance is negligible, thus 
Eq. ^ is complemented by continuous boundary condi- 
tions for Q at the ends of the wire. (The electronic energy 
distribution function for this problem has been discussed 
in Ref. H.) 

To accomplish this program technically we switch to 
the energy representation of Eq. Since the driving 
is periodic we can single out the i+ dependence for any 
operator!: 1(^1,^2) = E„^n(*-)e™"''*+- In the energy 
domain this implies that A{ei,e2) = J2n ^ni^i) 2nS{ei — 
£2 + nhfl), where A„(e — nhQ/2) is the Fourier transform 
of An{t^). It is thus convenient to define e = E + kMl, 
such that k in an integer, and —hn/2 < E < hil/2. 
One can then represents A in matrix form: AnmiE) ~ 
Am-n{E+nfl). In this way the time convolution between 
two operators becomes a simple matrix product. In the 
matrix representation the distribution function entering 
Eq. iQ for the right lead (z = L) takes a diagonal form: 
FR{E)nm = '^nmsign(i? + nil). For the left lead at zero 
temperature it reads instead: 



2l' 



1 — — OC 



Jii-A)J, 



rn—n — l 



(A) 



where A = eV° /Ml. In practical calculations it is enough 
to restrict the matrix size to \n\ < 3 A. 

We solve Eq. by representing the wire as a chain 
of TV — 1 ^ 1 chaotic cavities (where Q is uniform) with 
N identical barriers between them . This leads to 

a finite difference version of the Usadel equation: 



{Gk-i + Gk+i) + i 



N{N - l)eTh 



0, (14) 



where exh = h.D/L'^ is the Thouless energy of the 
wire. The operators e and have matrix represen- 
tations {e)n,n = E + nn and V''n,n{E) = ieV°{N - 
k){Sn,m-i - 5n,m+i)/{2N). Eq. ^ cau be solved by 
iteration |9j . Then photon-assisted shot noise can be ob- 
tained by evaluating the x-dependent current Ikix) = 
gNTi \^[Gk,Gk+i]fK} /8 at any one of the N barriers. 

Results for the differential noise dS/dVt versus AC fre- 
quency are shown in Fig. 13 We find that diffusion due 
to impurities induces on the photon-assisted noise a sim- 
ilar frequency dependence as the transmission through a 
chaotic cavity. Again a weak maximum is present with 
the main difference that the energy scale is set by the 
Thouless energy instead of the dwell time. 

In conclusion, we have shown that the frequency dis- 
persion of the photon assisted noise can be used to probe 
directly diffusion times in mesoscopic conductors. Our 
predictions can be verified by an experiment analogous to 
that described, for instance, in Ref. |3i where a chaotic 
cavity or a diffusive wire should be substituted to the 
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FIG. 3: Frequency dependence of the differential photon- 
noise for a diffusive wire for given values oi A = eV° /Ml. 
Here exh = TiD/L^ and the numerical calculation has been 
performed with 10 nodes. The magnitude of dS/dH is nor- 
malized to its value at H = 0. Curve (1) A = 1.0, (2) A = 2.0, 
(3) A = 3.0, (4) A = 4.0. 



quantum point contact. By carefully choosing the trans- 
parencies of the cavity or the length of the wire, one can 
match the Thouless energy {fi/T£i) with the range of fre- 
quencies that have been already investigated. A Thouless 
energy of 10 /^eV, that is typically realized in mesoscopic 
conductors, corresponds to Q./2tt = 2.4 GHz which is a 
frequency readily accessible in experiments. 
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